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We discuss scalar-tensor cosmology with an extra R~^ correction by the Noether Symmetry 
Approach. The existence of such a symmetry selects the forms of the coupling uj{ip), of the potential 
V{ip) and allows to obtain physically interesting exact cosmological solutions. 
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I. INTRODUCTION 
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CN ' Recent observational data indicate that ~ 70% of the today cosmological energy density is dominated by some form 

1—1 of "dark energy" which can be described, in the simplest way, by the cosmological constant A [ll 0, y, |j]- Such an 

, ^ , ingredient should explain the accelerated expansion of the observed universe, firstly deduced by luminosity distance 

measurements of SNela supernovae. However, even though the presence of a dark energy component is appealing in 

^N , order to fit observational results with theoretical predictions, its fundamental nature still remains an open question. 

■ Although several models describing the dark energy component have been proposed over the past few years, one of 

O , the first physical realizations of quintessence was a cosmological scalar field, which dynamically induces a repulsive 

<J • gravitational force, causing an accelerated expansion of the universe. 

Vh ' The existence of such a large proportion of dark energy in the universe presents a large number of theoretical 

OC problems. Firstly, why do we observe the universe at exactly the time in its history when the vacuum energy 

dominates over matter (this is known as the cosmic coincidence problem). The second issue, which can be thought 

^__l , of as a fine tuning problem,, arises from the fact that if the vacuum energy is constant, like in the pure cosmological 

J> ■ constant scenario, then at the beginning of the radiation era the energy density of the scalar field should have been 

t ' vanishingly small with respect to the radiation and matter components. This poses the problem that, in order to 

\l ■ explain the inflationary behavior of the early universe and the late time dark energy dominated regime, the vacuum 

y^ ' energy should evolve and should not simply be a constant. 

^^ . Some recent works have shown that the fine-tuning problem could be alleviated by selecting a subclass of 

C^^ ' quintessence models, which admit a tracking behavior ^], and in fact, to a large extent, the study of scalar field 
55 ■ quintessence cosmology is often limited to such a subset of solutions, fn scalar field quintessence, the existence condi- 
tion for a tracker solution provides a sort of selection rule for the potential V{(j)) (see [6| for a critical treatment of this 
question), which should somehow arise from a high energy physics mechanism (the so called model building problem). 
Also, adopting a phenomenological point of view, where the functional form of the potential V{(j)) can be determined 
from observational cosmological functions, for example the luminosity distance, we still cannot avoid a number of 
problems. For example, an attempt to reconstruct the potential from observational data (and also fitting the existing 
C^ ' data with a linear equation of state) shows that a violation of the weak energy condition (WEC) is not completely 
excluded 7], and this would imply a superquintessence regime, during which w^ < — 1 {phantom regime). However, 
it turns out that assuming a dark energy component with an arbitrary scalar field Lagrangian, the transition from 
regimes with w^ > — 1 to those with w^ < — 1 (i.e. crossing the so called phantom divide) are probably physically 
impossible since they are either described by a discrete set of trajectories in the phase space or are unstable ^, _9|. 
These shortcomings have been recently overcome by considering the unified phantom cosmology [lO| which, by taking 
a generalized scalar field kinetic sector into account, allows one to achieve models with natural transitions between 
inflation, dark matter, and dark energy regimes. Moreover, in recent works, a dark energy component has been 
modelled also in the framework of scalar-tensor theories of gravity, also called extended quintessence (see for instance 

P E M M MMMMMM)- 

It turns out that they are compatible with a peculiar equation of state w < —1, and provide a possible link to issues 
occurring in non-Newtonian gravity [12l |. In such theoretical backgrounds, the accelerated expansion of the universe 
is due to the effect of the non-standard form of the gravitational action. In extended quintessence cosmologies (EQ) 
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the scalar field is coupled to the Ricci scalar R in the Lagrangian density of the theory: the standard term ISttG* R 
is replaced by 167rF{(f>) R, where F{(j)) is a function of the scalar field, and G* is the bare gravitational constant, 
generally different from the Newtonian constant Gn measured in Cavendish-type experiments [ll| . 

Of course, the coupling is not arbitrary, but it is subjected to several constraints, mainly arising from the time 
variation of the constants of nature ^]. In EQ models, a scalar field has indeed a double role: it determines at any 
time the effective gravitational constant and contributes to the dark energy density, allowing some different features 
from the minimally coupled case [21]. Actually, while in the framework of the minimally coupled theory we have 
to deal with a fully relativistic component, which becomes homogeneous on scales smaller than the horizon, so that 
standard quintessence cannot cluster on such scales, in the context of non-minimally coupled quintessence theories 
the situation is different, and the scalar field density perturbations behave like the perturbations of the dominant 
component at any time. This is referred to as gravitational dragging ([l3|). 

On the other hand, the cosmic speed up can be simply explained considering some sub-dominant terms of geometric 
origin like i?^^, where R is the Ricci curvature scalar, which becomes dominant toward small curvature regimes (see 
e.g. [iH). In fact, it is possible to show that, by adding these terms to the Hilbert-Einstein action and varying with 
respect to the metric, such modified field equations naturally produce the observed cosmological acceleration. The 
simplest action of these models is: 




s = z^^ R-^W~9d'x (1) 



where Gn is the Newton's grav itational constant and /iq is a constant. The Palatini variation of this action is 
studied, for example, in Ref. [l9|. However, we need some additional ingredient to fit the observed data and physical 
constraints at every redshift so a modified scalar-tensor theory could be a more suitable candidate to achieve the 
whole observed dynamics. In this perspective, the investigation of theories like 
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^-7^)+^5'"'Va<0V.0-1/(0) 



-gd^x (2) 



is in order. Here (p denotes a real scalar field, non-minimally coupled to gravity while Lu{(f)) and V{4>) are the coupling 
function and a self-interacting potential, respectively. 

A physical criterion to achieve general cosmological solutions could be by the Noether Symmetry Approach which 
revealed a useful tool to fix the forms of the coupling and the potential [22], and, very recently, also the form of f{R) 
[23l . |24| . From a mathematical point of view, the method lies on the fact the presence of symmetries selects cyclic 
variables which allow to reduce the dynamics and then to integrate the equations of motion. From a physical point 
of view, any Noether symmetry is associated to some conserved quantity. This fact allows to select physically viable 
models (see for example [23]) and constitutes a criterion to select suitable effective Lagrangian (in particular, the 
forms of the coupling, of the self- interacting potential and the higher-order corrections). 

Specifically, in this letter, we work out the above action ([2]) searching for Noether Symmetries in order to see if the 
coupling, the self-interacting potential and the R^^ can be related in physically viable models. Besides, as we will 
see, this procedure allows to exactly integrate the equations of motion. 

The letter is organized as follows. In Sect. II, we search for Noether symmetries for the above action selecting 
the coupling and the potential. Sect. HI is devoted to find out the cosmological solutions and the discussion of the 
various sub-cases. Concluding remarks and conclusions are reported in Sect. IV. 

II. THE NOETHER SYMMETRY 

In order to search for Noether Symmetries, it is convenient to recast the action ([2]) by redefining (j) = (p'^ and 
Mo ~ ~Mj that is 



S 



f^ (i? + ^) + Mv)9'"'^^V^.V - V{v^)] ^^d\ (3) 



fa a? k\ 
Using the Friedmann-Robertson- Walker (FRW) metric, the scalar curvature takes the form R = &[ — \ — ^H — ^|, 

\a a^ a''' J 

where a{t) is the scale factor of the universe and the dot denotes the derivative with respect to time, with k = ±1,0 

. Now, using the Lagrange multipliers method [26[ , one can rewrite the action (3) as follows 



5* = 



^'{R+L^ + Va>(¥>) - Vi^) + Ai(i?-6- + - + - 



(4) 



where, the scalar curvature R and scale factor a are considered as two independent variables, and Ai is a Lagrange 
multiplier. The parameter Ai can be determined by varying the action with respect to R , that is 

Now, in order to apply the Noether symmetry approach, one can easily show that, in a FRW manifold, the Lagrangian 
related to the action (4) takes the point-like form [22|| 

C = 2a^(f'^iJ,q + &{p.q^ - l){2a'^ip(up + ip^a^a) + 12 y.ip^a^qaq - 6(f'^ kad^q'^ - 1) + a^(4w((^)(^^ - V{(f)) (6) 

where q = R^^ is a new variable. This means that we are considering an effective theory with two scalar fields. The 
corresponding Euler-Lagrange equations are given by 

ifiq^ - 1) \2^H^ + ^ + ^^) + 3 rf-"^ + 3'^(<^)<^ - 3 W^ + M^)^H + — — = (7) 



i^^q' 1) 



^2" ' 2a2 ' " ' dt 



2ip^H + ip'i-H^ + — +H+ ^^) 



+ ^V{v) - uj{^)^^ + M'/''^ + 2fi^qq{2^ + H^)^0 (8) 



with the constraint, derived from the definition of the scalar curvature R, 

6(2H' + H+^\=^ (9) 

Here iJ = — is the Hubble parameter. Eqs. (7) and (8) are equivalent to the second order Einstein equation and to 

a 
the Klein-Gordon equation, respectively. Finally, one can choose the initial conditions of these field equations such 

that the energy function associated with the Lagrangian (6) vanishes, that is 

^ .dC .dC .dC „ ^ , „, 

Ec^d— + q—+^—-C^Q, 10 

Oa oq dip 

or explicitly 

W - 1) (wff + \^^H^ + 1^) + \^{^)^^ + ^V{^) + fi^^Hqq - ^^lq^^ = (11) 

which corresponds to the {0, 0} Einstein equation. Now, let us introduce the lift vector field X [30| as an infinitesimal 
generator of the Noether symmetry in the tangent space TQ{a, a, Lp, (/?, g, q} related to the configuration space 
Q — {a, q, if} as follows 

d d d ■ d ■ d ■ d 

X = A^ + B^+C-^+A^+B^+C^ (12) 

oa 0(p oq oa oip oq 

where A, B and C are unknown functions of the variables a, <p and q. The existence of Noether symmetry for the 
dynamics implies that the vector field X is non-trivial and the Lie derivative of the Lagrangian, with respect to this 
vector field, vanishes 



Explicitly, this condition leads to the following differential equations 

Ga'^iiLp'^qA - 3a'^V(Lp)A - Ba^-— + 'ina^ipqB+ 

dip 

+ 2a^Lp^^iC + Qkip^A{l - fiq"^) + 12katpB{l - fiq^) - 12fikaip'^qC = (13) 

3Lu{tp)A + Ba—- + Sifiq^ - 1)¥^T- + 2u;{ip)a^— ^ (14) 

dip dip dip 

/ a A r)n\ f r)C\ 

(/V - 1) U'^ + 2Ba + 2aip^— + 2a^^— + 2^qap C + a—- = (15) 

\ oa oa J \ oa J 

, 1 ..f^ . r. dA ^dA dB\ 2 , , ^dB / ^ dC\ „ 

[m -^) [^VJ^ + ^a + ap—^p^—^pa—), -^ -uj(ip)a^ —^ ^qap>\2C + ip—\ =0 (16) 

/ BA Bn\ / r)A r)C\ 

{m^ -I) {p}-^+a—\ +2fiqAp + 2fiqaB + napiC + q—+q—j =0 (17) 

.dA 2 , ^ dB , 2 ,^dA ^ 



and 



Putting (19) into (18) implies 
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3fxqip^—+2Lo{ip)a—^0 (20) 

By choosing A = Aoa"p™, B = Bo{q)a''ip'' and substituting them into Eq. (20), we get 

„ / N 3 mAn n , , . 

Bo(g) = -TM -q^ + ki (21) 

4 luq 

i^{^) = ^0^"-^+' (22) 
where Ao^i^Q and fci are constant. By substituting this results into (14) we get 

too = m = 1 a-iid s = 2 (23) 
Taking into account Eqs. (21), (22) and (23), we get the solutions 

uj{p,) - 1 (24) 

A = Aoa"ip (25) 

B^{-^^iAoq^ + kl)a'"^p,'' (26) 

An important remark is in order at this point. In the case /x = 0, such solutions are ruled out, if ip is massless, 
by gravity tests on Solar System. This is not true for fi y^ 0. In this case, the previously mentioned tests strongly 



constrain the allowed masses for ip and therefore the parameters in the potential V{ip). For a detailed discussion on 
the effective scalar field mass constrained by Solar System tests see [23, [23, [231 ■ 
In view of these solutions, Eqs. (15), (16) and (17) read 



and 



7 2fci 

2^ aT 



{n^ - I)Aoa>2 



fci 



3/ig- 
9 



a'V M^o'Z + M'^V' [C + q 



dC 







(3 + n + ^)--,V 



These equations are satisfied if 
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- — - — -^q'^Aoa"(f'^ + aipnq ( 2C + ip^ 



dC 
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2fci 
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2fiqqip> I C 



dC 
da 



q' = ql = -G 

C = /oa"-V 



(27) 
(28) 



/3o Ao9o /^Q 2 7 o 4/ci 
/o = h — ;— 3/igo - 7 - 2n - -— 



G = 
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n An V 77, y 



Po 



A, 2n(n^ + n - 1) - 5 + ^(6n^ -in- 16) + %(4:n- 12) 
2^ n2 + n-3+^(n-2) 

/in 



where G, /o, 9o and /3o are constant. In conclusion, the Noether symmetry for the Lagrangian (6) exists and the vector 
field X is determined by (25), (26) and (27) and (28) while the functional form of u:{ip) is given by (24). 
It is straightforward to obtain a general self-interaction potential from Eq. (13) as 



where, we have used the definitions 



V{ip) = \p^ + k2p^ 



I. ^ foA ^ aA ' 12/i/o 



(29) 
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I-2A2 



Ai = ( 2go - 3g3 + |^ ) A^ 



(30) 



. 1 2 fci 
^' = r"' 3Ao 

The existence of the Noether symmetry means that there exists a constant of motion. In this case, the conserved 
quantity corresponding to the Noether symmetry can be obtained using the Cartan one-form associated with the 
Lagrangian (6), that is 

dC^ dC ^ dC ^ 
Oc = -rrrda + -rrrdp + -rrrdq 
oa oip oq 

By contracting 9c with X one obtains the following required constant of motion 

Fq = ixOc = Aoipa'\{fiq'^ - l){12a'^pp + I2p^aa -f 72fia^p'^qq)+ 



+ 12/oa"+VVa9 + a""V^(fci - tMo9o)(12(m9^ ~ l)a^dp + 8a^ujip)ip) 



(31) 



III. THE COSMOLOGICAL SOLUTION 

Starting from (9) and (27), it is straightforward to get the foUowing general sohition for the scale factor 



a{t) = J6fcgo + aiexp( ) + a2exp(— ==) 

V V^Qo v3<?o 



where ai and a2 are arbitrary constants. In special case, by choosing k = a2 — and ai — Oq, this solution takes 
the standard de Sitter form 

a{t) — flo exp(at) (32) 

where a = — ; Clearly this is a singularity free solution evolving as an hyperbolic cosine. It shows accelerated 

2V3(7o 
phases for t — > ±oo so both inflationary and dark energy behaviors are easily achieved. 

Some interesting sub-cases can be obtained considering the field potential (29). For ^2 = 0, it takes the form 

V{(p) = X(p^. In addition, one can use the constant of motion (31) and the scale factor (32) to find a solution for (p{t). 

To this purpose, we rewrite (31) as 

Fo = a"+2(Ao(M^ - l)(12(/3'<yi + 12^^a) + Ufo^^'fiaqo + ^^{h - ^fiAoq^){12{fiql - l)a^ + 8^)) (33) 

and then 

^(t) = ipo cM-^ot) (34) 



with 



and 



{n + 2)a fFa 

■t'o = t; , "^0 = 



Uq 
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uo = a^'+'(12Ao(wo - l)(12^o + 12a) + Ufot^aqo + (fci - -WoAo)(12a(Aigo - 1) + 8^o)) 

It must be stressed that these results have been obtained due to the existence of the Noether symmetry, and one can 
easily check that these solutions are consistent with the corresponding field equations. In this case (fc2 = 0), solutions 
(32) and (34) satisfy the Eqs. (7) and (8) which now assume the forms 

2 1 A 

2(pH'^{fiql - 1) + -(/5 - -fiqoip + 2ifH + -^ = 

{^lql - 1) (^2^^H + ^H\^ + ^^^ + iA^2 _ ^2 ^ 

respectively, with the following definitions of the constants 

2 _ 1 348n + lOSn^ + IGn^ + 473 
^°^^ (4n + 4n2 + 73)(13 + 2n) ^ ' 

583n + 327^2 + IGn^ - 4^^^ + 1508 
~ 27go(4n + 4n2 + 73) 

Clearly n is a free parameter depending on the constant of motion. 

Another interesting case is for A = and A2 = |. The self- interaction potential takes the form VJLp) = fc2'/3^. 
As it is well known, this potential is widely used in the discussion of vibrations of polyatomic molecules [25| and it 
is widely used in chaotic inflationary models also if, strictly speaking, the V{ip) ~ ip^ potential is the prototype of 
chaotic inflationary potentials [22, l3l|) IsJ • However, we have to say that quartic (p^ potentials are almost ruled out 
by current observations (see, for instance, |33jV 

In this case, the general solution of the equations of motion is achieved numerically while particular solutions have 
a power-law form as discussed in ,34] for non- minimally coupled theories without the R~^ correction. 



IV. CONCLUDING REMARKS 

We have explored the conditions for the existence of a Noethcr symmetry in a scalar-tensor theory of gravity, with 
an extra R"^ term, in which the coupling function and the generic potential are unknown. The motivation for this 
study is that we want to construct cosmological models capable of achieving inflationary and dark energy phases. To 
this goal, we need two fields leading the two eras which, in our case, are (p and q = R~^. 

We have shown that the existence of the symmetry fixes the coupling and the self-interacting potential which have 
physically interesting forms. Furthermore, it allows to achieve exact cosmological solutions which are singularity free 
and suitable to mimic inflationary and dark energy behaviors. 

A more physically appealing model should consider the role of standard perfect fluid matter and should fit also the 
dust dominated phase [23] , but in this cases the Noether symmetry cannot always achieved. 

However, also if we have considered a phenomenological model, the important lesson of this research is that, as 
shown also in other contexts [22, [2J| , the Noether symmetry is a powerful approach to select physically motivated 
solutions. 
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